EPJ manuscript No. 

(will be inserted by the editor) 



a^ : 

ON 

> " 
o 



C/3 



o 
o 



> 

ON 

O 



Critical properties of the metal-insulator transition 
in anisotropic systems 

Frank Milde**, Rudolf A. Romer, Michael Schreiber, and Ville Uski 
Institut fiir Physik, Technische Universitat, D-09107 Chemnitz, Germany 

Received: date / Revised version: date 

Abstract. We study the three-dimensional Anderson model of localization with anisotropic hopping, i.e., 
weakly coupled chains and weakly coupled planes. In our extensive numerical study we identify and char- 
acterize the metal-insulator transition by means of the transfer-matrix method. The values of the critical 
disorder Wc obtained are consistent with results of previous studies, including multifractal analysis of the 
wave functions and energy level statistics. Wc decreases from its isotropic value with a power law as a 
function of anisotropy. Using high accuracy data for large system sizes we estimate the critical exponent 
sls 1/ = 1.62 ± 0.07. This is in agreement with its value in the isotropic case and in other models of the 
orthogonal universality class. 

PACS. 71.30.-|-h Metal-insulator transitions and other electronic transitions - 72.15.Rn Localization effects 
(Anderson or weak localization) - 73.20.Dx Electron states in low-dimensional structures (superlattices, 
quantum well structures and multilayers) 
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We study numerically the problem of Anderson localiza- 
tion [|] in three-dimensional (3D) disordered systems with 
anisotropic hopping. Previous studies using the transfer- 
matrix method (TMM) multifractal analysis (MFA) 
] 1^ and recently energy-level statistics (ELS) |^,|^ showed 
, that an MIT exists even for very strong anisotropy. The 

■ values of the critical disorder Wc were found to decrease 
,, by a power law in the anisotropy, reaching zero only for 

■ the limiting ID or 2D cases. The main goal of the present 
' paper is to determine the critical exponent v of this second 
. order phase transition with high accuracy. It is generally 

assumed that v only depends on general symmetries, de- 
, scribed by the universality class, but not on microscopic 
• details of the sample Thus, anisotropic hopping should 
' not change i/. Recent highly accurate TMM studies report 
' z/ = 1.54 ± 0.08 |], I' = 1.58 ± 0.06 0, ly = 1.61 ± 0.07, 
and 1/ — 1.54 ± 0.03 ||ll[ for isotropic systems of the or- 
thogonal universality class. But for anisotropic systems of 
weakly coupled planes, = 1.3 ± 0.1 and v = 1.3 ± 0.3 
was found [||. We found in a recent high precision ELS 
study V = 1.45 ± 0.2 § for the same model. To clar- 
ify this situation, we compute the localization length by 
means of the TMM with high accuracy for large system 
sizes and apply a finite-size scaling (FSS) analysis which 
takes into account corrections to the usual one-parameter 
scaling ansatz |l^ . The resulting value of the critical ex- 
ponent 1/ = 1.62 ± 0.07 confirms the recent high accuracy 
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estimates. Thus the anisotropic Anderson model belongs 
to the same universality class as the isotropic model. 

Another interesting aspect of anisotropic hopping be- 
side the question of universality is the connection to exper- 
iments which use uniaxial stress to tune disordered Si:P 
or Si:B systems across the MIT |l2[|l|,|l|,^. While ap- 
plying stress, the distance between the atomic orbitals 
reduces, the electronic motion becomes alleviated, and 
the system changes from insulating to metallic. Thus, al- 
though the explicit dependence of hopping strength on 
stress is material specific and in general not known, it is 
reasonable to relate uniaxial stress in a disordered system 
to an anisotropic Anderson model with increased hopping 
between neighboring planes. 

In the experiments, a large variation of the value of the 
critical exponent v has been observed with suggested val- 
ues ranging from 0.5 |l^ over 1.0 1.3 Q, up to 1.6 
[ p^ . Possibly this "exponent puzzle" ||l^ is due to other 
effects in the experiments such as electron-electron in- 
teraction [|4| or sample inhomogeneities [^,|l^,|l^ which 
are usually ignored in the original formulation of Ander- 
son localization. Furthermore, the extrapolation of finite- 
temperature conductivity data down to temperature T = 
is open to debate and should perhaps be replaced [051 
|l8| by application of the dynamical scaling approach [p,9| . 
Another interesting question is, whether applying uniaxial 
stress is equivalent to changing the dopant concentration. 
We note that for non-universal properties such as the value 
of the conductivity, it was shown that stress and concen- 
tration tuning lead to different T dependencies close to 
the MIT ill. 
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2 The anisotropic Anderson model of 
localization 

We use the standard Anderson Hamiltonian 



H 



\m 



(1) 



i) weakly coupled planes with 

^ „ = 1 , = 1 - 7 
and (ii) weakly coupled chains with 



anisotropic transport 



1 



■ 7 



(2) 



(3) 



This defines the strength of the hopping anisotropy 7 G 
[0,1]. For 7 = we recover the isotropic case, 7=1 
corresponds to independent planes or chains. Note that 
uniaxial stress would be modeled by weakly coupled chains 
after renormalization of the hopping strengths such that 
the largest t is set to one in equation 



3 Transfer-matrix method in anisotropic 
systems 

We study the localization length A, describing the expo- 
nential decay of the wave function on long distances. We 
compute it using the TMM [|[||o[|ltl for quasi-lD bars of 
cross section M x M and length L ^ M. The stationary 
Schrodinger equation Htf' = is rewritten in a recursive 
form: 



{El - H,)/4 -1 







tf'i, Hj, and are wave function, Hamiltonian matrix, 
and transfer matrix of the ith slice, respectively. Unit and 
zero matrices are denoted by 1 and and tb is the hopping 
integral along the bar axis. We consider the band center 
E = 0. Given an initial condition (^^) equation (^) allows 
a recursive construction of the wave function in the bar 
geometry by adding more and more slices. A(M, W) is 
then obtained from the smallest Lyapunov exponent of 
the product TiTi_i • • • T2T1 of transfer matrices p^ , 
where the length L of the bar is increased until the desired 
accuracy of A is achieved. With increasing cross section of 
the bar the reduced locahzation length Am = A(M, W)/M 
decreases for localized states and increases for extended 
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with orthonormal states \i) corresponding to electrons lo- 
cated at sites i = {x,y,z) of a regular cubic lattice with 
periodic boundary conditions. The potential energies 
are independent random numbers drawn from the inter- 
val [—W/2,W/2]. The disorder strength W specifies the 
amplitude of the fluctuations of the potential energy. The 
hopping integrals tij are non-zero only for nearest neigh- 
bors and depend on the three spatial directions, thus tij 
can either be t^^ ty or t^. We study two possibilities of 
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Fig. 1. Reduced localization length for coupled planes with 
7 — 0.96 for parallel (top) and perpendicular (bottom) orienta- 
tion of the transfer-matrix bar. Accuracy is 1% (large symbols) 
or 0.2% (small filled symbols connected by lines to guide the 
eye). 



states. Thus it is possible to determine the critical disorder 
at which Am is constant from plots of Am versus M. 

For the anisotropic systems there are two possible ori- 
entations of the axis of the quasi-lD bar: parallel and 
perpendicular to the coupled planes or chains. The local- 
ization lengths in the perpendicular direction are smaller 
than in the parallel direction by a factor of about 1 — 7 for 
coupled planes and (1 — 7)^ for chains Nevertheless, 
the critical disorder Wc should not depend on the orien- 
tation of the bar For strong anisotropics 7 > 0.9 this 
is difficult to verify numerically, as can be seen for the 
case of weakly coupled planes in figure |l|. For the parallel 
orientation of the bar we find the usual critical behavior 
of Am as described above. We deduce a critical disorder 
Wc ~ 7 for this case. But for a perpendicular orientation 
of the bar the behavior of Am versus M is different as 
can be seen in the bottom part of figure 0. There are two 
striking features. First, Ai\i oscillates for small W and M 
between smaller values for odd and larger values for even 
M. Second, the characteristics of Am as function of M 
changes from localized (with positive slope) at small M to 
extended (with negative slope) at larger M for < 7. Let 
us consider for instance the data for W = 6. For M < 11 
Am decreases with M, which is typical for localized states. 
Up to M « 25 the data still decrease, but the slope tends 
to zero. For M > 25 it starts to increase, indicating ex- 
tended behavior. Therefore one has to extend the calcula- 
tion at least to M — 35 to find the correct critical disorder 
in this case. For smaller M, Wc would be systematically 
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underestimated even when applying the FSS procedure. 
We remark that, e.g., the computation of the data point 
for W = 8, system size M — 36, and accuracy of 0.2% 
takes several weeks on a 400MHz Pentium II machine. 
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Fig. 2. DOS for M weakly coupled planes of size M xM with 
W = 3 and 7 = 0.96 for an odd and an even M. 
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Fig. 3. Anisotropy dependence of Wc for coupled planes as 
computed by TMM in parallel (||) and perpendicular (_L) di- 
rection, previously by MFA and recently by ELS We 
also added a fit to TMM data of Ref. ||| (dashed line). Note 
the large systematic error explained in section ^ for the TMM 
in perpendicular direction at 7 = 0.99. 



We attribute these features of Am at least partially to 
a structured density of states (DOS) p{E) at these large 
7 and relatively small W. We show an example in figure 
g. The structure comes from very small (M x M) planes 
in the bar which are very weakly coupled in the perpen- 
dicular direction. The coupling between the planes is so 
small for 7 > 0.9, that p{E) is nearly equal to the DOS 
of an ensemble of uncoupled 2D systems g . In such small 
2D systems the relatively weak disorder is not sufficient 
to completely smear out the peaks in the DOS of the or- 
dered system. Thus, at E' = there is a peak for even but 
a dip for odd system sizes as can be seen in figure ^ In 
our opinion, for the TMM in perpendicular orientation, M 
has to be at least so large that all the finite size structure 
in p{E) has vanished in order to get reliable results. For 
the TMM in parallel orientation, smaller M are sufficient, 
since the planes or chains extend along the bar so that the 
DOS is smoothened. 



4 Computation of the critical properties at 
the MIT 

4.1 Anisotropy dependence of Wc 

Depending on the quality of our available data, we com- 
pute the critical disorder with different methods. The re- 
sults are shown in figures ^ and ^. Particularly for the per- 
pendicular orientation we estimate Wc from plots of Am 
versus M as in figure ^ As described above, a constant be- 
havior of Am for large system sizes indicates Wc- For data 
without the described features due to a structured DOS, 
i.e., in the parallel orientation, we plot the disorder depen- 
dence of Am for several system sizes as in figure |[ The 
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Fig. 4. Anisotropy dependence of Wc for coupled chains as 
computed by TMM in parallel (||) and perpendicular (_L) di- 
rection, previously by MFA and recently by ELS [0]. We 
also added a fit to TMM data of Ref. Q (dashed line). 



transition is indicated by a crossing point of the AMiW) 
curves. We use FSS for high quality data as described 
in the next subsections. Our results for Wc are in good 
agreement with results from ELS [|,@] and MFA g. The 
power-law dependence on anisotropy Wc = 16.5(1 — 7)^^ 
is confirmed. Using all data from MFA ELS, and the 
present TMM, we find /3 = 0.25 ±0.05 and (3 = 0.60 ±0.08 
for coupled planes and chains, respectively. The latter de- 
viates shghtly from the CPA result [|| /3 = 0.5. 
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4.2 Finite-size scaling 

The MIT in the Anderson model of localization is ex- 
pected to be a second-order phase transition [^|2^ . It is 
characterized by a divergent correlation length ^oo{W) — 
C\W — Wc\^'^ , where v is the critical exponent and C is 
a constant To construct the correlation length of the 
infinite system from finite size data Am pJ^,|20|,pT|, 
the one-parameter scaling hypothesis [E3 is employed, 



Am = /(AZ/Coo) 



(5) 



All A M are expected to collapse onto a single scaling curve 
/, when the system size is scaled by foo- In a system with 
MIT such a scaling curve consists of two branches cor- 
responding to the localized and the extended phase. One 
might determine v from fitting obtained by a FSS pro- 
cedure Q . But a higher accuracy can be achieved by fit- 
ting directly the raw data We use fit functions |0 
which include two kinds of corrections to scaling: (i) non- 
linearities of the disorder dependence of the scaling vari- 
able and (ii) an irrelevant scaling variable with exponent 
—y. Specifically, we fit 



Am = foiXrM'/n + M-yhiXrM^'n 



(6) 



1=0 



Xr{w) =W + Y^ hnW'' (7) 
n=2 



with w = {Wc — W)/Wc and a„i and expansion coef- 
ficients. Choosing the orders n-c and mr of the expansions 
larger than one, terms with higher order than linear in the 
W dependence appear. This allows to fit a wider W range 
around Wc than with the previously used linear fitting 
[^. The linear region is usually very small. The second 
term in equation ^ describes the systematic shift of the 
crossing point of the Am{W) curves visible, e.g., in 

figure H and its inset. This correction term vanishes for 
large system sizes, since the irrelevant exponent y > Q. 

For the nonlinear fit, we use the Levenberg-Marquardt 
method as in Ref. Q]. It minimizes the statistics, 

measuring the deviation between model and data under 
consideration of the error of the data points. We estimate 
the quality of the fit by the goodness of fit parameter Q 
p5| . It considers £^rid the number of data points and 
fit parameters. For reliable fits it should lie in the range 
0.01 < Q < 1 We check the confidence intervals ob- 
tained from the Levenberg-Marquardt routine by a Monte 
Carlo and a bootstrap method |2^. Additionally, we test 
whether the fitted values of Wc^ v, and y are compatible 
when using different expansions of the fit function, i.e., 
different orders and Wr 0. 



4.3 Determination of u 

We estimate v for coupled planes with strong anisotropy 
7 = 0.9, where we have the most accurate data. A parallel 
orientation of the transfer-matrix bar is used in order to 
avoid the problems discussed in section ^ Compared to 



< 




Fig. 5. Reduced localization length for coupled planes with 
7 = 0.9 and M = 5, 6, . . . , 17. The lines are fits of the data 
according to equations (^) and with rir = 3 and nir = 2. In 
the inset we enlarge the central region without the data points 
to show the shift of the crossing point. 



the perpendicular direction, the convergence of the TMM 
is much slower and the computing time to achieve a cer- 
tain accuracy increases remarkably. We computed Am up 
to M = 17 with 0.07% accuracy for 8.1 < < 9, for 
W = 8, 9.1, and 9.2 the accuracy is 0.14%. As we show in 
figure o we find a clear signature of an MIT, a crossover 
from increasing to decreasing behavior of Am with grow- 
ing M when disorder changes from 8 to 9.2. The lines for 
constant M do not cross exactly in a single point. In the 
inset, a small systematic shift is clearly visible. Thus, we 
include the second term of equation (0) when fitting the 
data. All fits reported in table |^ describe the data very 
well. This is expressed by the large values of Q > 0.7 and 
can also be seen in figure || where we show the data and 
the fit functions for an exemplary set of parameters n,-, TTir. 
The corresponding scaling function and scaling parameter 
are displayed in figure ^ and its inset. All data collapse 
almost perfectly onto a single curve with two branches. In 
connection with the divergent ^oo, this clearly indicates 
the MIT. Wc also tried to use smaller orders of the ex- 
pansions than in table ^ but then it was not possible to 
fit the data in the whole W interval with the desired high 
quality. 

When comparing the spreading of the fitted Wc and v 
values in table l| with their confidence intervals for the case 
that all system sizes are used in the fits (open symbols), 
the error estimates appear to be slightly too small. The 
95% confidence intervals of the smallest and largest Wc 
value do not overlap. We thus conclude Wc = 8.63 ± 0.02. 
In figure |^ we show the fitted v values and their confidence 
intervals together with results from our recent ELS study 
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Fig. 6. Scaling function and scaling parameter, shown in the 
inset, corresponding to the fit in figure |^. The symbols distin- 
guish different W values of the scaled data points. 

using the same fit method. The characters A to E de- 
note ELS results from different combinations of W and M 
intervals. Despite the high accuracy of the data of 0.2% to 
0.4% for large system sizes up to M — 50 and large Q val- 
ues, the results for scatter strongly. The 95% confidence 
interval apparently do not describe the correct error in 
this case. But for the TMM data, the error estimate for v 
seems to be appropriate. We emphasize the importance of 
having very accurate data for high system sizes as prereq- 
uisite to obtain rehable critical exponents. Furthermore, 
it is necessary to compare the results of different fits to 
get reasonable error estimates. 

In order to test for a possible systematic trend in the 
finite size behavior, we have repeated the fits neglecting 
the smallest system sizes M = 5,6. This is denoted by 
filled symbols in table 0. The results do not change, only 
the error increases. We summarize our result for the crit- 
ical exponent as = 1.62 ± 0.07. This is different from 
1/ = 1.3 ± 0.1 and larger than ly = 1.3 ± 0.3 obtained 
previously from data with an accuracy of about 2% for 
system sizes up to M = 15 and 17 for the parallel and 
perpendicular direction Since we use more accurate 
data with slightly larger system sizes we expect our result 
to be more reliable. Furthermore, u ^ 1.6 is in good agree- 
ment with high accuracy TMM studies for the isotropic 
case [p|, [LO|p^ . For comparison, we have added the results 
of Ref. |10(| to figure |7[ In our recent ELS study H we 
obtained v — 1.45 ± 0.2. As in other ELS studies p%p8t, 
the critical exponent is smaller than deduced from highly 
accurate TMM data. However, within the error bars, that 
result is consistent with our present finding. In Ref. [Q 
a trend towards larger v was found when the data from 
smaller samples were neglected. A further increase of v 
can be presumed if the system size could be increased fur- 
ther. We believe, that for large enough system sizes TMM 
and ELS will give the same results. 
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Fig. 7. v and its 95% confidence intervals from the fits of the 
Am data as reported in table ^. Large characters denote fits of 
ELS data from Ref. Q . The • data indicate results of Ref. Q 
for the isotropic case 7 = 0. Solid and dashed lines indicate 
the error bounds for the present result v — 1.62 ±0.07 and the 
ELS resuh = 1.45 ± 0.2, respectively. 

5 Summary 

We have studied the metal-insulator transition in the 3D 
Anderson model of localization with anisotropic hopping. 
We used TMM together with FSS analysis to characterize 
the MIT. Our results confirm the existence of an MIT for 
anisotropy 7 < 1 for weakly coupled planes and weakly 
coupled chains and the power law decay of the critical 
disorder with increasing anisotropy found in studies using 
TMM H, MFA g, and recently by ELS §,§. In these 
anisotropic systems there are two possible orientations of 
the transfer matrix bar. We have shown that the critical 
disorder Wc is, as expected, the same for both possibili- 
ties. But we remark that for strong anisotropy 7 very large 
system sizes are necessary for the perpendicular orienta- 
tion in order to find the correct Wc- This is in part due 
to the small size of the weakly coupled planes or chains in 
the bar which results in a structured DOS. 

For the case of weakly coupled planes with 7 = 0.9 
and parallel orientation we computed Am with 0.07% ac- 
curacy for system widths up to 17 x 17. Using a method 
to fit the data |l^ which considers corrections to scaling 
due to an irrelevant scaling variable and nonlinearities in 
the disorder dependence of the scaling variables we have 
deduced a critical exponent — 1.62±0.07. This is clearly 
larger than v = 1.3 ± 0.1 obtained previously Q for the 
same system. Since this result was obtained from less ac- 
curate data (« 2%) and slightly smaller system sizes, we 
believe that the previous error estimate is too small. Even 
from highly accurate ELS data (0.2% to 0.4%) and sys- 
tem sizes up to 50 the error estimate is twice as large: 
v = 1.45 ± 0.2 0. We have shown that large system sizes 
and high accuracies are necessary to determine the crit- 
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ical exponent reliably. Our result is in good agreement 
with other high accuracy TMM studies for the orthogonal 
universality class |^,10,ll|,^. These numerical estimates 
of V seem to converge towards v « 1.6. Experimentally it 
is of course even more difficult to determine the exponent 
V of the Anderson transition. Recent attempts of dynam- 
ical temperature scaling have shown that the statistical 
accuracy of the experimental data is less than in the nu- 
merical studies ||l^,^,|9l, but there also seems to be a 
trend towards larger values of v |l3,n5|. 



We thank K. Slevin and T. Ohtsuki for communication of their 
results prior to publication. This work was supported by the 
DFG within SFB 393. 



24. D. J. Thouless, Phys. Rep. 13, 93 (1974). 

25. W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. 
Vetterling, Numerical Recipes in FORTRAN, 2nd ed. 
(Cambridge University Press, Cambridge, 1992). 

26. K. Sle vin and T. Ohtsuk i, Phys. Rev. Lett. 78, 4083 
(1997), |2ond-mat/9704192 



27. E. Hofstetter , Phys. Rev. B 57, 12763 (1998), cond- 



mat/9611060 



28. I. K. Zhare keshev and B. Kra mer, Phys. Rev. Lett. 79, 
717 (1997), |cond-mat/9706255 . 

29. S. Bogdanovich, M. P. Sarachik, and R. N. Bhatt, Ann. 
Phys. (Leipzig) 8, 639 (1999). 



References 

1. p. W. Anderson, Phys. Rev. 109, 1492 (1958). 

2. Q. Li, C. Soukoulis, E. N. Economou, and G. Grest, Phys. 
Rev. B 40, 2825 (1989). 

3. L Zambetaki, Q. Li, E. N. Economou, and C. M. Soukoulis, 



Phys. Rev. Lett. 76, 3614 (1996), |cond-mat/9704107 , 

N. Panagiotides and S. Evangelou, Phys. Rev. B 49, 14122 

(1994). 

F. Milde, R. A. Romer, and M. Schreiber, Phys. Rev. B 
55, 9463 (1997). 

F. Milde and R. A. Romer, Ann. Phys. (Leipzig) 7, 452 
(1998). 

F. Milde, R. A. Romer, and M. Schreiber, submitted to 



Phys. Rev. B (1999), cond-mat/9909210 



8. B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 
(1993). 

9. A. MacKinnon, J. Phys.: Condens. Matter 6, 2511 (1994). 
10. K. Slevin and T. O htsuki, Phys. Rev. Lett. 82, 382 (1999), 



gond-mat/9812065 



11. P. Cain, R. A. Romer, and M. Schreiber^ Ann. Phys 



(Leipzig) 8, SL33 (1999), cond-mat/9908255 



12. M. A. Paalanen and G. A. Thomas, Helv. Phys. Acta 56, 
27 (1983). 

13. H. Stupp, M. Hornung, M. Lakner, O. Madel, and H. 
V. Lohneysen, Phys. Rev. Lett. 71, 2634 (1993). 

14. S. Bogdanovich, M. P. Sarachik, and R. N. Bhatt, Phys. 
Rev. Lett. 82, 137 (1999). 

15. S. Waf fenschmidt. C. Pfle iderer. and H. v. Lohneysen, 



(1999), |cond-mat/9905297 



16. T. F. Rosenbaum, G. A. Thomas, and M. A. Paalanen, 
Phys. Rev. Lett. 72, 2121 (1994). 

17. H. Stupp, M. Hornung, M. Lakner, O. Madel, and H. 
V. Lohneysen, Phys. Rev. Lett. 72, 2122 (1994). 

18. K. M. Itoh, M. Watanabe, Y. Outuka, and E. E. Haller, 
Ann. Phys. (Leipzig) 8, 631 (1999). 

19. D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 
(1994). 

20. J.-L. Richard and G. Sarma, J. Phys. C; Solid State Phys. 
14, L127 (1981). 

21. A. MacKinnon and B. Kramer, Phys. Rev. Lett. 47, 1546 
(1981). 

22. A. MacKinnon and B. Kramer, Z. Phys. B 53, 1 (1983). 

23. E. Abrahams, P. W. Anderson, D. C. Licciardello, and 
T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 



Frank Milde et al.: Critical properties of the metal-insulator transition in anisotropic systems 



7 





M 






X 


Q 


VVc 


u 


y 


A 


5 


■•17 


3 


1 


306.2 


0.789 


8.62 ±0.01 


1.65 ±0.04 


1.56 ±0.27 


< 


5 


■ ■ 17 


2 


3 


309.8 


0.745 


8.64 ±0.01 


1.59 ±0.04 


1.31 ±0.23 


V 


5 


■ ■ 17 


3 


2 


303.0 


0.815 


8.63 ±0.01 


1.64 ±0.04 


1.51 ±0.27 


l> 


5 


■ ■ 17 


3 


3 


300.7 


0.829 


8.63 ±0.01 


1.64 ±0.04 


1.55 ±0.27 


▲ 


7 


■•17 


3 


1 


218.6 


0.995 


8.64 ± 0.03 


1.66 ±0.07 


1.34 ±0.77 




7 


■ • 17 


2 


3 


211.7 


0.998 


8.65 ±0.02 


1.60 ± 0.05 


1.34 ±0.47 


T 


7 


■ ■ 17 


3 


2 


209.2 


0.999 


8.64 ± 0.02 


1.62 ±0.07 


1.38 ±0.51 


► 


7 


■ ■ 17 


3 


3 


208.9 


0.998 


8.65 ±0.03 


1.59 ±0.12 


1.24 ±0.58 



Table 1. Fit parameters and estimates for Wc and f with 95% confidence intervals from fitting Am for coupled planes with 
7 — 0.9. The symbol in the first column is used in figure m 



